MULTIPLE ANALOGUES OF BINOMIAL COEFFICIENTS AND 
RELATED FAMILIES OF SPECIAL NUMBERS 
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Abstract. We construct multiple ijrt-binomial coefficients and related multi- 
ple analogues of several celebrated families of special numbers in this paper. 
These multidimensional generalizations include the first and the second kind of 
gt-Stirling numbers, gt-Bell numbers, qt-BernouUi numbers, gf-Catalan num- 
bers and the gt-Fibonacci numbers. In the course of developing main proper- 
ties of these extensions, we prove results that are significant in their own rights 
such as certain probability measures on the set of integer partitions. 



1. Introduction 

Many distinct sequences of special numbers are investigated and their properties 
are explored in number theory. Various generalizations of such families of numbers 
have also been studied extensively for most cases. An important class of generaliza- 
tions of special numbers is their one dimensional, one parameter q-extensions. In 
this paper, we give multidimensional basic qi-generalizations for several collections 
of such numbers including the binomial coefficients and the classes of numbers as- 
sociated with the names of Stirling, Bernoulli, Catalan, Bell and Fibonacci. We 
also point out how to construct multiple ordinary a-generalizations for the same 
number sequences. The definitions of these numbers and the properties they satisfy 
show great variety which makes this research area very interesting. 

Among many mathematicians who contributed to this line of research, L. Carlitz 
appears to be the first to study the g-extensions for several families of special 
numbers given in this paper. Many have made significant contributions since then 
investigating properties of ^-generalizations, their applications and connections with 
other types of numbers. We will give references to some of these successful efforts 
and important work in section |4] below. 

The present paper takes a step in generalizing the one dimensional q-special 
numbers to multiple qi-special analogues. These generalizations are given in terms 
of the qt-hinomial coefficient defined in section [3] in the form 

'A\ glMlt2n(M) + (l-n)|/^| r (qp-')^ 



^^Jq,t (9^" l<z<j<n ^ 't'^-t'^ 

where A and fj. are n-part partitions and g, i G C. The W^^ function that enters 
into the definition is a limiting case of the well-poised BCn Macdonald function 
Wx. We first give a brief review of this remarkable family of functions and its 
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one parameter generalization uj\ called well-poised BCn Jackson coefficients. The 
symmetric rational functions W\ and uj\ are first introduced in the author's Ph.D. 
thesis |15j in the basic (trigonometric) case, and later in [16' in the more general 
elliptic form. 



2. Background 

We start with the definition of the g-Pochhammer symbol (a; q)a, where q, a 6 C, 
which can be defined formally by 

(2.1) (a;q)„:= 



(ag";9)c 



in terms of the infinite product (a;g)oo '■— Y^iLo{^ ~ '^'t)- Note that when a — m 
is a positive integer, the definition reduces to (a; q)„i = Y\jk=ai^ ^ o.q'^)- 
An elliptic analogue is defined [THl [3H] in the form 



m — 1 



(2.2) {a;q,p),n:= Ueiaq"') 



k=0 



where a S C, m is a positive integer and the normalized elliptic function 0{x) is 
given by 

(2.3) 0{x) = e{x;p) := {x;p)ao{p/x;p)oo 

for x,p E C with \p\ < 1. The definition is extended to negative m by setting 
{a;q,p)m — ^/{aq"^',q,p)-m- Note that when p — 0, {a;q,p)m reduces to the 
standard (trigonometric) ij'-Pochhammer symbol. 

For any partition A = (Ai, . . . , A„) and < G C, define [41] 

n 

(2.4) {a)x = (a;g,P,t)A l[{at^'';q,p)x,- 

fc=i 

Note that when A = (Ai) = Ai is a single part partition, then [a]q,p,t)\ — 
{a;q,p)xi = (a)Ai- The following notation will also be used. 

(2.5) (ai,...,afc)A = {ai, . . . ,ak;q,p,t)x := {ai)x . . . {ak)\- 

Now let A = (Ai, . . . , A„) and /i — (/xi, . . . , fin) be partitions of at most n parts 
for a positive integer n such that the skew partition A//i is a horizontal strip; i.e. 
> Ml > A2 > /i2 > ■ • • A„ > A'n > Ki+i — Mn+i — 0- Following [T7], we define 

(2.6) Hx/^{q,p,t,b) 



(qM.-W-i + ltJ-»-l)^^._^_;,^.(5A.+A, + 1^2-j-»5)^^_^_ 



n 

l<i<j<ri 

(gA.-p,_l+4j--.-l)^^_^_^^ I (g^'+^^+4^-^-^&)^,_,-A, 
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and 



, , (x ^ .ax)\(qbx/t,qh/(axt)),, 

(2.7) w,,,ix, q,p, t, a, b) := H,^,iq,p, t, b) • ^-i^ 



■^1 9{bt^-^^) {bqt-^^),,+,^^, ■ 

where q,p, t, x,a,b £ C. The function Wa/^(j/, zi, . . . ,zt, q,p, t, a, b) is extended to 
^ + 1 variables y, zi, . . . , G C through the following recursion formula 

(2.8) Wx/^,{y, zi, Z2, . . . , zf, q,p, t, a, b) 

= Wx/^iyt^'^;q,p,t,at'^\bt^) W^/^{zi,.. . ,ze,;q,p,t,a,b). 

We will also need the elliptic Jackson coefficients below. Let A and /i be again 
partitions of at most n-parts such that A//i is a skew partition. Then the Jackson 
coefficients wa/^j are defined by 

{x-\ax)x {qbr-^x,qb/axr)^ 

(2.9) WA/^(a;;r,9,p,t;a,6) — — ^ — ^ 

[qbx,qb/ax)x (x ,ax)f^ 

{rM~'t^'n, A f e{br-H^-'^q^^^) 
(g6r-2,gt"-i)^ Ai\ 0{br-H^-^') > 

■Wf,{qH'^"^;q,p,t,bt^-^",br-h^-") 

where x, r, q, p, t,a,b (1 C. 

Note that W\/fj,{x;q,p,t,a,b) vanishes unless A//i is a horizontal strip, whereas 
uix/fj,{x;r;a,b) = iux/ ^(x;r,q,p,t;a,b) is defined even when A//i is not a horizontal 
strip. 

The operator characterization [17 of i^Jx/n yields a recursion formula for Jackson 
coefficients in the form 

(2.10) ujx/r{y,z;r;a,b) Yujx/t,{r^''y;r;ar^'' ,br'') uj^/^{z;r;a,b) 

where y = {xi,. . .,Xn-k) S C""'' and z = {x„-k+i,. . . ,a;„) S C''. 

A key result used in the development of the multiple special numbers, the cocycle 
identity for ujx/fj,, is written in |17| in the form 

(2.11) uj^/^{{uv)~^-,uv, q,p, t; a{uv)^, buv) 

= X! ^i'/A(^'~^^^:'7,P,^;a(ww)^,&wu)wA/p(M"^w,g,p,t;aM2,&M) 

/iC ACi/ 

where the summation index A runs over partitions. 

Using the recurrence relation (j2.10l) the definition of ujx/f^{x]r]a,b) can be ex- 
tended from the single variable x € C case to the multivariable function llIx/i_i{z; r; a, 6) 
with arbitrary number of variables z = {xi, . . . , Xn) S C". That LUx/^iz', r; a, b) is 
symmetric is also proved in 1171 using a remarkable elliptic BCn io'P9 transforma- 
tion identity. 
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Let the Z-space V denote the space of infinite lower-triangular matrices whose 
entries are rational functions over the field F = C{q,p, t, r, a, b) as in [17]. The con- 
dition that a matrix u ^ V is lower triangular with respect to the partial inclusion 
ordering C defined by 

(2.12) ^J- <^ X ^J.i < Xi, Vi > 1. 
can be stated in the form 

(2.13) — 0, when fi % X. 

The multiplication operation in V is defined by the relation 

(2.14) (uw)am ^ uxuVv^_, 

for u^v £V . 

2.1. Limiting Cases. The limiting cases of the basic (the p = case of the elliptic) 
W functions Wx/ ^{x\q,t,a,h) = W\/fj^{x;q,0,t,a,b) will be used in computations 
in what follows. To simplify the exposition, some more notation will be helpful. 
We set 

(2.15) W^}^{x;q,t,s) \iiR Wx/^{x;q,t,a,as) 
and 

(2.16) WlJ^{x;q,t) \\ms\^\~\'^\W^lix;q,t,s) 
and finally, 

(2.17) W^j^ix; q, t) := lim W^)^{x- q, t, s) 

The existence of these limits can be seen from (p = case of) the definition (|2.7p , 
the recursion formula (12.81) and the limit rule 



(2.18) lim al'^l(x/a)„ ^ (-1)1^1 xl^li-"^^)^"^'^') 

a— >0 

where = X^Li M« and n(/i) = ~ and n{^') = YJl^i C^')- 

These functions are closely related to the Macdonald polynomials [32j |36] and 

BCn abelian functions [37] . 

We now make these definitions more precise. Let H\/^(q, t, b) — Hx/fj,{q, 0, t, b), 

and define 

(2.19) Hx/^,{q,t):^rmiHx/^,{q,t,b) 



n 

l<i<j <n 



By setting 6 = as in the definition of Wx/^ function, and sending a — ?■ we define 



the family of symmetric rational functions W?'/ in the form 



(2.20) W^lix;q,t,s) := \imWx/f,{x;q,t,a,as) 



,-n(X) + \t,\+n(t,) rr ( ^)\iqs/{xt))^ 
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for X e C. Using (|2.8p we get the following recurrence formula for W^^^ function 
(2.21) W^liy, z; q, t,s) = J2 ^A/.(2/^"'; 1, t, st-') ^(^z- t, s) 

u-<X 

where y G C and z £ as before. 

^A/p and Wl)^ 



Similarly, the W^^t and W^^, are defined as follows. 



(2.22) Wl)^^{x;q,t):^ s\^\~\^\W^)^{x-q,t,s) 



The recurrence formula for M^ft function turns out to be 

\/ n 

(2.23) M/0^(y, z- q,t)^Y. t'^^^^-^'^^W^^^M^': q, t) Wl]^{z- q, t) 

for y € C and z G C^. In the same way, we define 

(2.24) W;^j^{x; q, t) lim W^)^(x; q, t, s) 

_ ,-n(A) + |M|+n(M) , C„ f] ^^ 

The recurrence formula for Wf] function becomes 

A/^ 



(2.25) Wlj^iy, z: q,t)^Yl KjM^'' 9' K/^^^■' 

i/^A 

where again y E C and z G 



3. gt-BiNOMiAL Coefficients 

A common property for all types of special numbers we develop in this paper is 
that they are closely connected with binomial coefficients. Therefore, we start with 
the definition of g<-binomial coefficients which will be proved using a multiple ana- 
logue of the terminating version of gi-binomial theorem. We first derive a multiple 
terminating 24'i sum as a limit of Jackson's 8(^)7 summation formula. 



Theorem 3.1. For an n-part partition X, we have 

(3-1) ^^-E'^'^'^'^'^^r^ n (rl^Y^ 

•M^;''(g^i''(");g,t,s~ir~i) 

where s, x,t,q E C. 

Proof. The one variable basic (i.e., p = 0) version of the Jackson sum [17] 
(3.2) ujx{s~'^x;s;as'^,bs) = U!x/f_,{s~'^;s;as'^,bs)ujf^{x;s;a,b) 

^CA 
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can be written explicitly in the form 

^2 g-j {sx~'^,asx)x _ is,as)\ {bt^~",qb/as)^ 

{qbx,qb/ax)\ ^^{qb,qb/a)\ (qt"'~\as)^ 



Substituting b = da and sending a, d — > gives the desired result. □ 

The terminating qt-binomial theorem follows as a corollary. 
Corollary 3.2. For an n-part partition X, we have 

(3.4) {x)x = J2 (-l)l^lgl^l+"(^')r('^'+(i-")l''l 



l<'i<j<n 

where x,t, q € C 

Proof. Replace x — )■ sx~^ and s — > in identity p.ip in that order, and send 
s — using p.l8|) to get the identity to be proved. □ 

This latter identity is a multiple analogue of Cauchy's q-binomial theorem |24) . 
Using p.4p we give the definition of a multiple analogue of the binomial coefficient 
as promised. 

Definition 3.3. Let A and /i be n-part partitions. Then the gt-binomial coefficient 
is defined by 

where q, t G C. 

Note that with this definition we can write the terminating gt-binomial theo- 
rem (13. 4p in the form 

(3.6) {x)x^^{~l)^^^q''^'''h-"^f'^(^^ xl'^l 

Note also that setting t — q'^ and sending q — )■ 1 yields multiple ordinary a- 
binomial coefficients. Below we extend this definition further to be valid not only 
for partitions A and /i, but also for any 7i-tuples of complex numbers A G C" and 
/J, G C". For n = 1, the definition reduces to that of the one dimensional g-binomial 
coefficients 

[q)n-k{q)k 

which are also known as the Gaussian polynomials that are studied extensively in 
the literature including but not limited to the works in[il[51[Ml|51[^|551[^[Tl]. 
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Some of the main properties of the binomial coefRcients readily generalize to the 
multiple case. For example, the identities 

(3.8) '""-tit) 0-E(-l)'(fc) 

fc=0 ^ ^ k=0 ^ ^ 

have the following multiple analogues. 
Theorem 3.4. For an n-part partition X 

(3.9) (-I)a - (Z"'^'^t""^^^ 
and 

(3.10) = (-l)l''lg"(^')t-"('') 

Proof. These identities follow immediately from p.6p by setting z = — 1 and z = 0, 
respectively. □ 



Before investigating multiple analogues of other binomial identities, we first 
note that the definition p.3p makes sense even for generalized partitions fi — 
(MIj /^2, • ■ • , /^n) such that /ii > /i2 > • • • > /in with possibly negative parts 
starting with some i £ [n] := {1, 2, . . . , n}. 



Theorem 3.5. The W\ function is well-defined defined for the generalized parti- 
tions A. In fact, the evaluation of Wx function that enters into the definition of 
qt-binomial coefficient can be computed for n-tuples of complex numbers X S C" 
and /i € C". The limiting W^'' , W^^ and W^^ functions are also defined in the 
general case. 

Proof. Let A be an n-part partition with A„ ^ and < /c < A„ for some integer 
k, and let z = {zi,. . . , Zn) G C". We have shown [16j that 

(3.11) W,^{z;q,t,a,h) 

ilbP-'-''hk l\ {qbzM<lb/iaz,))k'^'-'"^''^ '"'^ '^-^ 



Among other applications, this identity extends the definition of Wx function to 
general partitions. For a generalized partition A, the index A — fc" on the right hand 
side will be a standard partition. 

The duality formula for Wx functions from [16] may be used to compute the 
special evaluation of Wx that occur in many application in this paper even for 
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A G C". The formula can be stated as 



(3.12) Wx {k-^q^t^; q, t, k'^a, kb) 



n 

l<z<j<n 



(fc)A(fcat"-i)A 

(iJ-'+i)A,-A,(5a't2"-*-J)A,+A, 

{qbt^-^),iqb/a% 



{h-^qh^;q,t, h'^a',hb) 

n 

l<'i<j<n 



where /c = a't" and /i = at" ^/6. Substitute these latter relations between 
parameters into the identity and factor out one term 



corresponding to the dominant weight Ai from the W function on the left hand 

side, and one term 



for ui from the W on the right hand side to get 



(3 15) 9-^-"(fc)M(g/fc).^.(fca*^^"-^^g-0A. ^, .^_r ^ 



{qbt^-^)x{qb/a)x ^ f (f^-*)A.-A, (g/c6r-*-^)A,+A, 



(fc)A,(fct-l)A'(fcat«-l)A ^J/.^^ I (t^-+l)A.-A,(5fc6t"-^-^ + l)A.+A, 

{{bt^-''/a)qH^; q, t, a^kf^/b, a^-^) 



{(at^-^/b)q-^\akP(^-^)q^^U /„^A.^. „ . „2,,„-l 



{qbt^'^),{qt^-^/k), -j-j f (*■'-')..-., M2"-'-i-l),^+,^ 



where W denotes the W function after the dominant factors arc taken out and A' 
denotes the partition (A2, A3, .... A„). Setting k = 1, moving entries to the right 
hand side, and multiplying back both sides by the dominant factors for the left 
hand side gives 

(3.16) Wxiq^t^-^q^t.a^b) 

q^''''{q-''')x^ it-')yiat^-')x rr ( {P''+')x^-x,{qbt^-'-^+^)x^+x, 



n 



(q)^, (qfet"-l)A(<zV«)A I (t^'-OA.-A,(gfot"-^-^)A.+A, 

^^"'""(bi'-V+Ar'"'^^" i(bt'-Va)qH^^-hq,t,aH--'/b,at--') 



{qbt^-%{qt^-^), -pr ( iP-%^^,^{qat^"-'-J-^),^ 



{at^-^bUat^-^), ^^^^^{{p-i+^U..,iqat^^-i-dU+.^ 
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This formula allows us to extend the evaluation of the Wx function on the left hand 
side from an n-part partition A to any A G C", because A appears as a variable 
on the right hand side. The front factors that involve A can be extended to the 
complex case by the definition of g-Pochhammer symbol (|2.1[) . 

Similar evaluations for the limiting cases Wf and W^^ can easily be com- 

puted from this result. □ 

We now turn to the study of main properties of gt-binomial coefficients. 

Theorem 3.6. For an n-part partition A and generalized n-part partition fi, we 
have 

(3.17, f!) ^ (t) - 1. 



-.K, V0"/„,, 

where 0" is the n-part partitions whose parts are all zero. In addition, 
(3.18) (^) =0 



when % \ or (f^ pL. 



Proof. Recall that the normalizing coefficients for Wx function is written [16] in 
the form 

(3.19) Wx{q^t'''-^^;q.t,a,h) 

^ A [ {qbt^-\qt--'^)xM^''-^%X. .(n+l-2k)X, 

H I ((a/6)^"-^ar-'=)A,(g6^"+l-2'=)2A, 

•^"/^'^'^^ n iqt^-^)x.^xiat^+--^-^)x.^X- 

Setting b ~ as and sending a ^ and s — >■ cxd (after multiplying by s'^') gives 



(3.20) W^f(gV(");g,i) = (<?i"-^)At("-^)l''-^"(')g-l'l J] r 

l<i<j<n ^'^ A, -A, 

Substituting this evaluation into (^)^ ^ gives 1 as desired. 

It is easy to see from the definition of Wx function that Wp(x; q, t, a,b) = 1 when 
u = 0". Therefore W^„ (x; q,t) = l and (^) = I. 

Finally, the fundamental vanishing property of W functions states (TB] that 

(3.21) W^{qh^;q,t,a,b) = 

when ^ 2 A. Hence W^(g^i'^("); g, t) = and therefore (^)^ = in that case. Note 
also that if < then (^) ~ due to the (gi"^^)^ factor in the denominator. □ 

We introduce some notation before writing a multiple analogue of the important 
recurrence relation 

^n + 1\ / n\ f n \ 

\<k< 



k J \kj \k~l 

for classical binomial coefficients. Let — (0, . . . , 1, . . . , 0) be the n-dimensional 
vector whose ?-th coordinate is 1 and all others are 0. Let A* = A + whenever A* 
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is a partition. Finally, we write A h fc to denote that A is a partition of the positive 
integer k. 

Theorem 3.7. For an n-part partition A and k < \X\, we have 



(3.22) ^g"(^')t 



where /i C A' and t,i> C X. 

Proof. The proof is a simple application of the gi-binoniial theorem (I3.6p . Since 

= (1 - xt^-''q^'){x)x, we get 

(3.23) ^ (-i)lMlg"(p')i-«(f)/'^''j 



(1 - xt^-'g^O ^(-l)l^lg"(^')i-"(^ 



by a double application of p. 61) . This may be written as 

I A' I 



(3.24) ^^(-l)l^lg"('^')t-"('^)('^'') 



= (1 - xii-'g^O E E(-l)''''^"^''^^""^'^ (i) 
Note that the coefhcient of x*' on the left hand side is 

(3.25) 

The coefficient of x^ that comes from the two pieces on the right hand side becomes 
(3.26) Y.^-lt\q^^-^'h''<-^(^\ + E (-l)l''l+lg"('^')+^'i-"('^)+l-'('^' 

Note also that /i C A* and t, C A, for otherwise the binomial coefficients vanish 
by Theorem (|3.6p above. Canceling out common factors gives the result. □ 



A multiple analogue of the symmetry property 



n\ n 



, I 1 , < A: < n 

kj \n — k J 

for classical binomial coefficients now follows. 

Corollary 3.8. For an n-part partition A and k < \X\, we have 
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Proof. This result follows from the gt-binomial theorem p.6|) through a similar 
argument used in the proof of the previous theorem. □ 

We will now write an analogue of the identity 

(3-28) (r^ = ^Q)^"' 1^1 <^ 

^ ^ n=k ^ ^ 

Let A be an n-part partition and consider the C(g, t)-space Pa of all polynomials 
of degree less then or equal to |A| in a single complex variable x. Note that both 
sets /3i — {x^ : k = 1, . . . , |A|} and — {(a^)^ : ^ A} form bases for Px- In 
fact, qf-binomial theorem stipulates just this, providing a change of basis formula 
between the two bases. 

Similarly, we view the terminating 2(/>i sum p.ip as a matrix representation of the 
shift operator acting on /32 , and prove an simpler version of the cocycle identity |16j 
for oja functions as follows. 

Theorem 3.9. For n-part partitions v and fi, we have 

(3.29) (sr), W^;^((z''i*("); g, t, isr)-H"~') 

■MAVl^;''(gV(");g,t,r-4"-i) 

Proof. Write the identity p.ip in the form 

(3.30) (..)a=E'?'^'^^"^'^'7^ n {(SSr^ 

^CA 't' l<^<J<n ^ ^'^^ >t''-f'^ 

•W^;'('zV(");g,t,s-ir-i)(x)^ 

Starting with the basis /32 , apply the shift operator by a factor of r followed by a 
shift by s. This double shift can be achieved by a composite shift by a factor of sr. 
Writing this argument explicitly using the identity above and simplifying gives the 
desired result. □ 

Using this result, we now write a multiple analogue of the binomial identity (|3.28p 
given above. Here (a)oo" denotes n"=i("^^~*)oo- 

Theorem 3.10. For 7i-part partitions v and fi, we have 

where max{|gzi(^*^"~-'^) | : i £[n]} <1. 
Proof. First send r — > in the (I3.29P to get 
(3.32) iy;(g''t^(");g,t) 

= y ^|A|i2„(A) (f). TT I (^^!2)^^Vab( .^.H ^ l^„-l) 
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Next we send — > cx). Note that setting v — A:", where fc" — {k,k, . . . , k) denotes 
the n-part partition whose parts ah equal k, simphfies the W function evaluations 
in the identity. To make this precise, we recall the Weyl denominator formula |16j 
which states that 

(3.33) W^{xt^^^^;q,t,a,b) 



Set b = cs^iii-" and a = ci^^^™ ([sjaj, and send c -> to get 

Multiplying both sides by (s~^t"~^)l''l , and sending s ^ further implies that 
(3.35) W;^(2:i*(");g,0 = (-l)l''lxl'^lr('')q-l^|-"(^')(:E-i)^ ]J *^')a'.-m 



l<?'<j <n 



On setting x — q^ ^ sending fc — > oo in the last two evaluations and substituting 
them into the identity (I3.32p . we get 

iiP'^+'),.-,, ^ ^ t2„(A)^(„-l)|A|,|A| 
^ '°° l<i<j<n ^ ADp ^ 

_ TT r (g^^'^)A.-A, (^'"^')^^-A4^st(A^^(n). ^) 

Replacing s by gives the identity to be proved upon verifying the convergence. 

The fact that the resulting infinite series converges follows from a multiple ana- 
logue of the dominated convergence theorem introduced in [17] . Consider the mul- 
tiple series 

(3.37) 

where is the sub-alcove {A G Z" : fc > Ai > . . . A„ > 0}. The theorem states 
that if the pointwise limit h\ := limfe_>oo h\{k) exists for all A G L^, and we can 
find for each A such that \hx{k)\ < m!{ for aU fc > Ai, and that the series 
'Yl,\eL+ "^A convergent, then the original series converges. 

That the pointwise limit exists on both sides in (|3.32l) is already verified above. 
Note that the index of the W function inside the summand is a fixed parti- 
tion. Note also that, except the powers of g, t and s, all other factors in the 
summand of (I3.32p can be put into the form {uq°')aol {vq°')oo using the definition of 
g-Pochhammer symbol (j2.ip . Standard theorems on infinite products and sequences 
imply that such factors are bounded when a is a non-negative integer, m, w G C, 
and that v is such that the denominator never vanishes. This is because 



(3.38) lim 



a— )-oo 



(ug")c 



1 
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when |g| < 1. Therefore it foUows that for a constant Ch that may depend only on 
q,t and s, and is independent of k and A, we've found 



(3.39) 



m'l = Ch 



f2n(A)-(n-l)|A|^|A| 



Ch 



i=l 



(2i-ri-l)Ai^Ai 



FinaUy, we need to show that X]a6L+ '^a convergent which may be verified using 
the multiple series ratio test. Let be defined as above. Then, we see that 



(3.40) 



A+e; 



TO, 



(2i-n-l) 



qz 



for each i G [n] where possible (i.e., when A + is a partition). Therefore, 
Sagl+ "^a converges when 



max{ 

ie[n] 



(2i-n-l) 



qz 



}<1. 



The original series also converges under the same condition. 
Next, we write a multiple analogue of the identity 



(3.41) 



E 

k—m 



for gi-binomial coefficients. 

Theorem 3.11. For n-part partitions v and fi, we have 



(3.42) 



t 



-n(M)„n(M')Izllii ("^ - 

(-1). w,.. 



= E ^ 

f-tC.XGv 



Proof. The cocycle identity ^16) for w functions 

(3.43) WA/^((sr)"\ sr, as^, hs) = ^ a;A/,.(s"\ s, as^, hs) u^/^,{r~^-,r, a, b) 

may be written explicitly in the form 

{rs)x(asr-^)x {qb)x{qb/a)x (qb/as)^ {ar^^)t. 



(3.44) 



{qbr ^)x{qbr/a)x {s)x{as)x {asr i)^ {qb/a)f^ 

(qb/as), {bt^~'% {r),{ar-^)u 



E' 



n 



(1 - 6t2-2z) W 



(as)^ ((?i"-i)^ ((j6r-i)^(gfer/a)^ 

(gt^-*).,-., {bt''-'-^),^+. 



n 

l<^<j<n 



□ 
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Setting b — azrs/q, and sending a — >■ and r — >■ gives 
(3.45) ifik (fksl^lT^^^((z''t'(");q,i) 



We will use a multiple analogue of Bailey's lo^g transformation formula from |16j 
to transform the series on the right hand side. The lo^g transformation can be 
written explicitly in the form 

(3 46) {v)Aa's^'"-^)u {qbr)^{qbr/a')^ {Av^s^)^ {qbrs/Av)^ 



E 



{qbrs/v)u{qbrv/sa')^ {s)i,{a's)^ {qbrs"^ /Av^)^ (Av/sr)^ 
{qbrs"^ /Av^)x {{brsv^^t^^^)x [r s / v) \{Av / sr) x {qbr/a's)\ {a'sv~^)\ 
^ {Av/s)x (gt"-i)A {qb)x{qbr^syAv^)x {a's^v-^)x {qbr/a')x 

fr\ {l~brsv-H^-^\^^^) x^\ yr \ (gt^--)A.-A, {brsv-H^-^-^)x^^x^ 

11 \ (1 _ 5^,^-1^2-2^) W ) ] ^VV^y[^v-^-^)x^_^^{brsv-H'^-^-i)x,^x^ 

■ Wx{t, brst^-^", brsv-H^-"; q^t^^""^ ; q) VK^(t, brsy-H^-^'' ,bf-''; g^i**"^ g) 
{r)x{Avyrs^)x {qbrs/Av)x {qbr/a's)x {brt^~'')x {a'sv-^)x 



^ {qb)x{qbr^s^/Av^)x {Av/s)x {a's)x {qt''-^)x {qbrv/sa')x 

/j,C ACi/ 

Af (l-&rt^-V^O ^.,,2.-2^Aa TT f (gt^-)A.-A, (^r^^^^^^" ) 

(l_6rt2-2,) W ; /^^ll^J(gt. — 1),,_,^. (6rt2— .),,+,^. 

• Wxit, brst^-'^'^, brt^-'';q''t^^"'^;q) VK^(i, fert^^^n^ ^^l-n. q^t^{n).q) 
Sending 6 — > and setting A — Asr/v gives 



(3.47) |:^(_l)lMl^lMli-«(f)gn(M') 

(^)m 



A ^ l<«<i<^^ 

/j-C ACi/ 



(„l)|AU|A|^|A|^-n(A)gn(A') ^^^^^^ ^^^^^^^ 



/.iC ACzy 



We now apply the v = s and A = — s^^ case of the last transformation identity on 
the sum side of p.45p with z = — s^^. The result follows after some manipulations 
of terms and simplifications. □ 



MULTIPLE BASIC SPECIAL NUMBERS 



15 



A multiple analogue of the identity 
(3.48) 2"-in = X:fcQ 

k=l ^ ^ 

follows from previous theorem as a special case. 
Corollary 3.12. For n-part partitions v and fi, we have 

(-1). (^{^-q''')t'-'\ 



(3.49) t-"Mqni^')}i^ 



eiCACi/ \i=l ^ / ^ ^9.* 

Proof. This follows immediately by setting /x = ei in p.42p and using the special 
evaluations p.60p studied below. □ 

One of the most important families of discrete probability distributions is the 
binomial distribution whose density is given by 

(3.50) /(z;n,p):= (^"^p'(l-p)""' 

where < p < 1 and i = 0,1,..., n. Here the density f{i;n,p) denotes the 
probability that the "event i" occurs for fixed parameters n and p. Changing the 
roles of i and n — i would give an equivalent definition. 

We now define an analogous probability measure on the set of all n-part parti- 
tions contained in A under the partial inclusion ordering. 

Definition 3.13. Let A be an n-part partition. For any partition /i C A, the 
qi-binomial density function is defined by 

(3.51) g{p;X,z):= zl^l-l''l(z),, 
where p G R with < p < 1 . An alternative definition would be 

Jul 



(3.52) 



/(/i;A,z) ■.^t-'ni,)^2ni,')fX\ 

VJq,t 



It needs to be verified that the density function is non-negative for any /i C A, 
and that the total probability adds up to 1 when summed over all /i C A. This is 
what we verify next. 

Theorem 3.14. The multiple discrete density functions g{pL]\,z) and f(fi;X,z) 
are valid densities. 

Proof. Set a — qb/r in (|3.3I) and send 6—^0 and s to get 

(3.53) = y iiMi(-l)IH n(M') (X\ 

Setting z — \/x now and sending r — > cxd gives 

(3.54) i=E-'"''^'w.n 
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which shows that g{^; A, z) add up to 1. Similarly, shifting x to xr in p.53|) and 
sending r ^ instead gives 

(3.55) 1 = E zl'^lr2«(A')^2"(M') j"^' 

verifying total probability for the alternative definition /(/i; A, z). 

It is obvious from the definitions (|2.4p and (13.31) that both densities are non- 
negative when < q,t, z < 1 and z < t. □ 

For a fixed n-part partition v, let 5'j/ be the set of all partitions that are under 
v with respect to the partial inclusion ordering. That is, 

S"^ := {A : A C i/} 

It then follows from Theorem 13. 141 that the distribution 

(3.56) F{\;i.,z):=Y,z^'^^-^'^Hz)J'' 

defines a probability measure on S^,. An alternative distribution may be defined 
using the / density function defined above as well. The relation between this 
measure and the one given in [TO] on the set of all partitions is to be investigated 
in another publication. 

We will, however, define a multiple analogue of another important family of 
density function for the Poisson distribution in this section. It will be defined as 
a limiting case of the binomial as in the classical case. First we give two multiple 
analogues of the exponential function e^. 

Theorem 3.15. The following functions are multiple analogues of the exponential 
function . 



(3.57) E,iz) := {-zU^ = ^ 



n 

l<z<j<n 



an 



d 



I 2:^fi\f'^n{fj.) + {l-n)\fj.\ 

where max{|zi'^^*~"~-'^^ | : i £ [n]} < 1. 

Proof. Setting A — fc" in gt-binomial theorem p.6p . sending — )■ oo by using the 
identity (j3.34p and applying the multiple analogue of the dominated convergence 
theorem employed in the proof of the previous theorem, and replacing x by —z 
gives the expression for Eq{z). 

Similarly, sending a; — oo in p.ip . setting A = /c" as above and sending k ^ oo 
using the identity p.35p . and replacing s by z gives eq{z). Note in the latter case 
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that the convergence theorem requires the condition max 

to be satisfied. □ 
We point out only two obvious properties of these functions. Namely, that 

e,{z)E,{-z)^l 

and that 

lim lim Eq{z{l — q)) ~ , and lim lim en{z{l — q)) — 

Both of these properties follow immediately from their definitions given in the 
previous theorem. Other properties Eq{z) and eq{z) satisfy will be investigated in 
a future publication. 

We now give a multiple analogue of the gt-Poisson distribution. 

Theorem 3.16. The function 



(3.59) fi^i■,z)■.= Eq{-z) 



n 

l<i<j<n 



defines a valid density on the set of all partitions of length at most n. 

Proof. Similar to the classical case, we set A = fc" and send fc — > cx) in the qt- 
binomial density f{fj,; A, z) using the identity (j3.35p to get the gt-Poisson density 
/(/i; z) function as desired. 

That the total probability X^/^eP /(/^' ^) ^'^'^ '^P ^^'^ that /(/i; z) is always 
non-negative follow from the construction. □ 

Finally, we study certain special evaluations of the qt-binomials. It turns out 
that when either A or is a rectangular partition /c" = (fc, fc, . . . , fc), then (^) ^ has 

a closed form product representation. 

Theorem 3.17. For n-part partitions A and ^i, the special cases 

have closed form product representations. 

Proof. The A = /c" case of identity p. lip becomes 

^QR1^ w f + u\ TT (9^^^"^')2fc A {z^^)k{az,)k 

3.61 Wk^{z-q,t,a,h) = \\ ,,,-_i_2»\ \\ (h\(hii vT 

i<»<i<n ('^^^ hklJi{qhzi)k{qb/{aZi))k 

Set b — as and send a — s- to get 

(3.62) w^i!^{z-q,t,s)^\{}^^ 
Furthermore, multiplying both sides by s"*^ and sending s — oo gives 

n 

(3.63) W^l{z; q, t) - n(9'''^.)fc 
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after flipping certain factors using 

(3.64) (a)„ = (gi-"/a)„(7('^')(-ar 

Therefore we get 



k 



It is clear from the definition that 



(3.66) WrlJx;q,t)^l 



and 

(3.67) W:l^4x;q,t)=q-\l-x) 

for X ^ C. It then follows from the recurrence relation (|2.23p that 

n / ^ \ 

(3.68) w:J{z;q,t) = J2 - ^7't"-l = q-' [J2x, + f-' 

1=1 \i=l / 

for z = {xi, . . . ,Xn) S C". Multiplying by the front factors gives 
1 - g" . 

where [nL := is the so-called q- number or g-bracket. 

1 - g 

The proof of the last result follows immediately from the identity p.35p . Setting 
A = fc" and manipulating factors gives 

K \ _ ,2»(m) + (1-«)ImI TT ^ ^A'. 



(3.70) ^ ^2„(p) + (l-„)|,.| TT , 



n 

l<i<j <n 

as claimed. □ 



f (g^-'-')M.-p, (t^"-'+^)M,-P: 



In the light of last theorem, we give a definition for gt- number [z]qt extending 
that of a g-number defined above. The definition will be used in the next section 
in the discussion of the qi-Stirling numbers. First, we write an extension of our 
qi-binomial coefficients. 

Definition 3.18. Let z = {xi, . . . , Xn) € C" and fi be n-part partitions. Then the 
extended qt-binomial coefficient is defined by 

/^\ „|Ai|+2n(^i) + (l-n)|^| ( (nP~^\ 1 

where q,t ^ C It should be noted that this definition makes sense even for fi G C" 
by the virtue of Theorem 13.51 when z is an n-part partition. 

With this extension, a definition for gri-number may be written as follows. 
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Definition 3.19. Let z = (si, . . . ,a;„) e C". Then 

We also define the partition shifted generalization [z;/i]^t = of the qt- 

bracket in the form of a normalized W^^ function as follows. 



(3.73) [z;A^],t:=<zl''l+"(^')n 



i=l 



1 



(1 - qt"-^)f'i 

n { (^Z!r" ]w;;iqH^^-\q,t) 



l<i<j<n 



The evaluation (I3.72p in the definition follows from p.63p . Note also that when 
z = g^t'^f") for a single variable x G C, the /x-shifted qt-mimheT [z; fi]qt can be 
written as 

" r 1 

(3.74) [x; - <r^' ^ (9^; I/'?, 1/0. n I (1 _ 
by using the following application of p.64p 

(3.75) x\^'\{x-\q,t)^ = (-l)l^l(?"('')i-"('')(x;g-\t-i)^ 
to simplify the M^"^ function. Similarly, using p. 631) we get 

[q q I )k 



(3.76) [z;fc%t:=g"(^)n|- 



(1 - qt"-')'^ 



In particular, for fc = 1, we recover the definition p.72p above. 

Note also that (x; l/g, l/i)^„ with the reciprocals of q and t, corresponds to a 
multiple basic qi- analogue of the falling factorial Xn-= x{x — 1) ■ ■ ■ {x — (n — 1)). 

4. Multiple basic and ordinary special numbers 

In this section we give multiple basic and multiple ordinary analogues (or the 
qt- and a-analogues, respectively) of several celebrated families of special numbers 
including Stirling numbers, Bernoulli numbers. Bell numbers, Fibonacci numbers 
and pentagonal numbers. The definition and some immediate properties of these 
numbers are studied in each case. Their fascinating deep properties such as the 
recurrence relations they satisfy, closed form evaluations in certain special cases 
and their combinatorial interpretations are carried out in other publications. 

The multiple special number sequences we introduce in this section are indexed 
by partitions. For each class of multiple basic gt-special number, the corresponding 
multiple ordinary a-special number analogues can be found by setting t = q"' and 
sending q — 1. 

4.1. qt-Stirling numbers. The Stirling numbers of the first kind are defined to 
be the coefficients of the falling factorial Xn in the expansion 



(4.1) Xn=n\r\ =^si{n,k)x'' 



k=0 
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The g-analogue of these numbers are defined in |7] and their properties are studies 
in [35l [27l [30l |43] and others. We give the definition of the muhiple gi-Stirhng 
numbers generaUzing the one dimensional g-analogues as follows. 

Definition 4.1. For an n-part partition A, the gi-Stirling numbers of the first kind 
si(A, ^i) are defined by 



(4.2) [x;A],i = ^si(A,M)n[^i""tt 

fiCX i—l 

where x G C. Consider the infinite dimensional lower triangular matrix m whose 
entries are si(A,/^). That is, m\f^ = Si(A,/x). The gt-Stirling numbers of the second 
kind S2{X, /i) are defined to be the entries of the inverse m^^ of the matrix m. That 

is, S2{\,fi) := m^^. 

We now write an explicit formula for the gt-Stirling numbers generalizing one 
dimensional g-analogues as follows. The proof of the theorem makes clear that the 
definition makes sense by presenting an inverse for the matrix m. 

Theorem 4.2. For n-part partitions v and /i, an explicit formula for the qt-Stirling 
numbers of first and second kind Si(A,/^) and S2(A,/i) are given by 

n(iy')Y--2n(^) + (ri-l)|M| s, > , / \ 

(4-3) -1 - nr-i(i-g^"-)--^' ■ ^ ^''^^ 



(4.4) S2{iy,fi) ^ S2{i^,fi,q,t) 

-n(^')^2ri(iy) + (l-n)|i/| , \ / m , 

^= n" (1 - at-'V-^-A'^ • ^ (li-"(-,A,l/g,lA))t("-^)l^lx'(A,A.,g,0 
where u{X, fi) and v(X, ii) are defined by 



(4.6) v{X,n,q,t) 



l<'i<_7<?i 



MULTIPLE BASIC SPECIAL NUMBERS 21 



Proof. Wc write the M^-Jackson sum in the form 
(4.7) Wxix;q,p,t,at-^'',bt-'') 

iqbs-H-^)x{qbt-^s/a)x ■ I (iJ-*)A.-A,(g&t-^--')A.+A, 

Cn^"-iv.r««-if-"-iv. Ill n -^«-lfl-2^^ ^-^^ 



n 

l<i<j<n 



(tJ--')^.-^^-(g*^-')^.-^^. (6s-ig*-'-^)^.+^^.(b,s-it--^^-+2)^.+^^. 
(qt^— i)^,_^,(i^-^+i)^._^, (6s-it— ^■+i)^.+^,(g6s-it— ^+1)^.+^^. 



Set 6 = ar in this identity and send a — >■ 0, r — >■ (X) using the hmit rule p.lSp . and 
set z = xt^^""^ and simphfy using p.35p to get 



(4.8) (x;g-i,ri)A = ^u(A,M)a: 



where u(A,/i) is defined as above. Similarly, replace x by x/s in (14. 7p and send s 
to to get 

(4.9) = 5]^>(A,M)(x;g-i,t-i)^ 

/iCA 

where w(A,/i) is defined as in the theorem. Then it is clear, by a change of basis 
argument, that 

(4.10) J,A= ^^(^>A)z;(A,m) 

Note that the left hand side (|4.9p does not depend on q or Also, 

(4.11) [xr-i ; 1, t] = n i-^-t.' - n rii'A^^;.' 

i=l i=l ^ ^ ' 

for X G C. By sending g — > 1 in identity (|4.9p and substituting it in (|4.8p after 
replacing x by xt"~^, we get 

(4.12) (a;;g-i,t-i)A = ^(a;r-i;l,t-i)^ 



^ w(^., A) ( Ihn v{X, fi, 1/q, lA)) 



This gives the explicit expression for the gt-Stirling numbers of the first kind after 
multiplying both sides by appropriate factors coming from the definition p.74p of 
[x;X,q,t\. Similarly, using the inversion relation ()4.10p gives the explicit formula 
for the qi-Stirling numbers of the second kind, S2{v,lJ-)- O 

These identities reduces in one dimensional case (i.e., n = 1) to the formulas for 
the g-Stirling numbers given in [30] . 
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An obvious property that fohows from the inversion relation (|4.10p in the last 
proof is that, considered as matrix entries, the gi-StirUng numbers of the first and 
second kind form inverse matrices. That is, 

(4.13) (5yA = ^ si(iy, A)s2(A,^) = ^ 32(1^, X) si{X, 

Leaving the investigation of other properties to a future pubhcation, we point out 
another immediate property in the next theorem. 

Theorem 4.3. For an n-part partition X 

(4.14) si(A,A) = .S2(A,A) = 1. 

Proof. Setting b = ast"~^ and sending a 0, and then s — > in p.l9|) gives 

(4.15) W^f (QV(");q,0 



= (_l)|A|^-«(A)^-|A|-„(V)(^i«-l)^. -Q 

l<i<j<n ^'^ 

Substitute this in u{X, A) to get 

(4.16) u(A,A) = (-1)1^1 
Similarly, substituting (|3.20p into v{X, A) gives 

(4.17) v{X, A) = (-l)l^lg"(^') t-"(^) 
In addition, the front factors 

and its reciprocal reduce to )iT^n{\)^{n-i)\x\ ^j-^g^ — ^. This is precisely 

what we want to complete the proof, because the factors (/^"(a ) jjj ^(^\^ ^j^j 
v{X,X) becomes 1 in the limit as g 1. Finally, the intermediate term and the 
flips q 1/q and t l/t cancel with the factor i-2n(A)+(«-i)|A| ^ □ 

It should be pointed out again that setting t — q" and sending q ^ 1 gives the 
multiple ordinary a- Stirling numbers of the first and second kind. 

4.2. gt-Bernoulli numbers. Another very important family of special numbers 
is the classical Bernoulli numbers Bk that may be defined in terms of a recurrence 
relation. Namely, Bq = 1 and Bk, for all integers fc > 1, satisfy 

n-l / s 
fe=0 ^ / 

where n > 2. Generalizing this approach in terms of multiple basic qi-binomial 
coefficients, we now give the definition of gt-BernouUi numbers. The vector ei and 
the partition A^ are defined as above. 

Definition 4.4. For an n-part partition A, the 5<-Bernoulli numbers are defined 
by the recurrence relation 

(4.19) = H i""'^'?"^^'' (^^) /3m 

mCA ^ ^ 
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with the mitial condition that /?0" — /3(o,o,....o) = 1- 
Solving the equation (|4.19|) for (3x gives 

(4.20) /3a = _r(A)q-n(A') j"^'"! ' ^ i-n(M)^«(p') f^'\ 

which can be used to compute (5^ for ah partitions /i. 

In the one dimensional (i.e., n = 1) case, the definition of qt-Bernoulli numbers 
reduces essentially to g-BernouUi numbers defined in and studied by [311 [30l [11] 
[To] and others. In the limit as g — ?► 1 in the one dimensional case, the gt-Bernoulli 
numbers reduce to the classical Bernoulli numbers. 

We conclude this section by the following conjecture. 

Theorem 4.5. In the limit as t = q°' and q ^ 1, the qt-Bernoulli numbers /?> 
vanish for all n-part partitions A whose weight is |A| = 2k + 1 for any integer k > 1. 

A proof of this conjecture among many other deep properties of these numbers 
will be investigated in a future publication. 

4.3. qi-Bell numbers. Another important class of special numbers intimately re- 
lated to the binomial coefficients is the Bell numbers defined by 

n 

B„ := ^S2{n,k) 

k=0 

This definition also generalize readily to the multiple case as follows. 

Definition 4.6. Let A be an n-part partition. The qt-BeW numbers are defined by 
the relation 

(4.21) := J2 i""^''^9"^'''^S2(A, fi) 

where S2(A, /i) is the qt-Stirling numbers of the second type. 
The classical Bell numbers satisfy the recurrence 

Bn+l ^^[f^jBk 
k=0 ^ ^ 

Among other properties, a multiple analogue of this result for B\ in terms of ^ 

will be given in another publication. 

For n = 1 this definition reduces to g-Bell numbers are defined and studied 
in [31], [ig and others. 

4.4. g^-Catalan numbers. A gt-generalization of the Catalan numbers already 
defined in [55] and developed in [321 US and others. These numbers are indexed by 
positive integers (the weights of partitions) and reduce to the g-Catalan numbers 
under the specialization of the t parameter. The multiple gt-Catalan numbers 
defined in this paper appear to be related to the definition given [22], however 
the precise relationship is still to be given. The gt-Catalan numbers of this paper 
reduces, when n = 1, to the one-dimensional q-Catalan numbers that are defined 
in [5] and studied in [1 HSDl HO] , and others. 
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The classical Catalan numbers are defined by the relation 

(4.22, 

n + 1 \ n 

We now give the following multiple qi- analogue of these numbers. 

Definition 4.7. Let A be an n-part partition. Then the gt-Catalan number Cx is 
defined by 



q,t 



where = A + ei = (Ai + 1, A2, . . . , A„) as before. 

In the special case when X — is a, rectangular partition, we get 

^ ■ ) L.fe" ■ _ qk+Hn-l) LL (1 _ qktn-^) LI {qt^-r)^ 

using the evaluation (j3.65p from above. In the particular case when fc = 1, this 
reduces to 

(4-25) ^--riTr^^ 

which, as an empty product, becomes 1 when n = 1. 
Similarly, using p.69p we see that 



(4.26) C, 



ei 



2ei 



Other properties of gt-Catalan numbers will be investigated in a future publication. 

4.5. gt-Fibonacci numbers. Recall that the classical Fibonacci numbers F„ have 
the representation 

(4.27) F„«. Y. {t 

n—k+m 

In a similar way, we extend the definition to the multiple basic case as follows. 

Definition 4.8. For an n-part partition A, the gi-Fibonacci numbers are defined 
by the relation 



(4.28) Fai q2«(A.')i2(n-l 



)[ph2n(p)/'^ 



where the sum runs over all n-part partitions C A and C A such that their 
sum equals A. Here, the sum v -\- ^ oi two partitions v and fi is defined to be the 
coordinate-wise summation of the two. That is, 

which itself is a partition. 
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In one dimensional case n — 1, the qi-Fibonacci numbers are equivalent to the 
q-Fibonacci numbers defined in [39] and studied in [8], [1], [21], [12], [13], [28] and 
others. In particular, in the limiting case g ~> 1 for n = 1, both definitions generate 
the classical Fibonacci sequence. 

Among a myriad of fascinating properties of the Fibonacci numbers, the funda- 
mental recurrence relation 

Fn — Fn-l + Fn-2, 

which is used to define the Fibonacci sequence with the initial conditions -Fq = and 
Fi = 1, can be generalized to qt-Fibonacci numbers Fx using Theorem 13.71 This 
fundamental result and other properties will be given in a separate publication. 

We close this section by pointing out that the relation between the so-called finite 
form of the classical Rogers-Ramanujan identities and the g-Fibonacci sequence 
holds true for our multiple gt-Fibonacci sequence. More precisely, it is known [1] [12] 
that the q-Fibonacci sequence may be generated by 

(4.29) F„+^{x)^ ^'l'^-'''' 



k=0 



k 



in the limit as 9 — > 1 for x — I. The finite form of Rogers-Ramanujan identities is 
essentially a limiting case of the Watson transformation. The multiple analogue of 
the Watson transformation given in [17] can be written in the form 

iqb,qb/p2(T2)x qM^2n(t,) i'^2. P2,qb/ pi(Ti)^, 
{qb/a2,qb/p2)x {qt^-^,qb/(7i,qb/pi)f, 

• n lifZlt'" I W;:\qH^^-y,q,t,P2a2e-^/qb) 
l<i<j<n ^ ^'^^ ^^'-^^ J 

-^Mpia2P2) ^ ' l\\ (l-6t--) / 



{bP '\cr2,P2,cri,Pi)ti 



{qf^ ^,qb/(7l,qb/pl,qb/a2,qb/p2)^, 
The limiting case of this identity, as cri, pi, (72, P2 00, becomes 



(4.31) iqb)x J2 gl^l+2"('^')t-2"(M)^,lMl I 



l<2<j<n 



which is a multiple analogue of the finite (terminating) version of the multiple 
Rogers-Ramanujan identities studied in [T7]. The summand on the left hand side 
becomes identical to that in our definition of gt-Fibonacci numbers on making the 
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substititution b = q ^t^^" except that the only partitions v and ji that enter the 
sum in the definition are the ones that add up to A. 

5. Conclusion 

The classical sequences of special numbers have fascinated mathematicians for 
centuries, and great many properties of these numbers have been studied. Several 
one or two parameter extensions of these numbers have been developed as well. 
However, these extensions were mostly one dimensional, and multiple analogues 
of these numbers have been missing. This paper takes an important step and 
presents multiple (basic and ordinary) analogues of several important classes of 
special numbers and develops some of their basic properties. Needless to say, the 
greatly rich properties of these numbers can not be fully covered in a single article. 
In addition, no attempt is made to include other families of related special numbers 
in this paper. These ideas will be pursued in future publications. 
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